Transcriptomic dysregulation and autistic-like behaviors in Kmt2c haploinsufficient mice rescued by an LSD1 inhibitor.

Recent studies have consistently demonstrated that the regulation of chromatin and gene transcription plays a pivotal role in the pathogenesis of neurodevelopmental disorders. Among many genes involved in these pathways, KMT2C, encoding one of the six known histone H3 lysine 4 (H3K4) methyltransferases in humans and rodents, was identified as a gene whose heterozygous loss-of-function variants are causally associated with autism spectrum disorder (ASD) and the Kleefstra syndrome phenotypic spectrum. However, little is known about how KMT2C haploinsufficiency causes neurodevelopmental deficits and how these conditions can be treated. To address this, we developed and analyzed genetically engineered mice with a heterozygous frameshift mutation of Kmt2c (Kmt2c+/fs mice) as a disease model with high etiological validity. In a series of behavioral analyses, the mutant mice exhibit autistic-like behaviors such as impairments in sociality, flexibility, and working memory, demonstrating their face validity as an ASD model. To investigate the molecular basis of the observed abnormalities, we performed a transcriptomic analysis of their bulk adult brains and found that ASD risk genes were specifically enriched in the upregulated differentially expressed genes (DEGs), whereas KMT2C peaks detected by ChIP-seq were significantly co-localized with the downregulated genes, suggesting an important role of putative indirect effects of Kmt2c haploinsufficiency. We further performed single-cell RNA sequencing of newborn mouse brains to obtain cell type-resolved insights at an earlier stage. By integrating findings from ASD exome sequencing, genome-wide association, and postmortem brain studies to characterize DEGs in each cell cluster, we found strong ASD-associated transcriptomic changes in radial glia and immature neurons with no obvious bias toward upregulated or downregulated DEGs. On the other hand, there was no significant gross change in the cellular composition. Lastly, we explored potential therapeutic agents and demonstrate that vafidemstat, a lysine-specific histone demethylase 1 (LSD1) inhibitor that was effective in other models of neuropsychiatric/neurodevelopmental disorders, ameliorates impairments in sociality but not working memory in adult Kmt2c+/fs mice. Intriguingly, the administration of vafidemstat was shown to alter the vast majority of DEGs in the direction to normalize the transcriptomic abnormalities in the mutant mice (94.3 and 82.5% of the significant upregulated and downregulated DEGs, respectively, P < 2.2 × 10-16, binomial test), which could be the molecular mechanism underlying the behavioral rescuing. In summary, our study expands the repertoire of ASD models with high etiological and face validity, elucidates the cell-type resolved molecular alterations due to Kmt2c haploinsufficiency, and demonstrates the efficacy of an LSD1 inhibitor that might be generalizable to multiple categories of psychiatric disorders along with a better understanding of its presumed mechanisms of action.

Chromosomal DNA sequences of the Pacific saury genome: versatile resources for fishery science and comparative biology.

Pacific saury (Cololabis saira) is a commercially important small pelagic fish species in Asia. In this study, we conducted the first-ever whole genome sequencing of this species, with single molecule, real-time (SMRT) sequencing technology. The obtained high-fidelity (HiFi) long-read sequence data, which amount to ~30-folds of its haploid genome size that was measured with quantitative PCR (1.17 Gb), were assembled into contigs. Scaffolding with Hi-C reads yielded a whole genome assembly containing 24 chromosome-scale sequences, with a scaffold N50 length of 47.7 Mb. Screening of repetitive elements including telomeric repeats was performed to characterize possible factors that need to be resolved towards 'telomere-to-telomere' sequencing. The larger genome size than in medaka, a close relative in Beloniformes, is at least partly explained by larger repetitive element quantity, which is reflected in more abundant tRNAs, in the Pacific saury genome. Protein-coding regions were predicted using transcriptome data, which resulted in 22,274 components. Retrieval of Pacific saury homologs of aquaporin (AQP) genes known from other teleost fishes validated high completeness and continuity of the genome assembly. These resources are available at https://treethinkers.nig.ac.jp/saira/ and will assist various molecular-level studies in fishery science and comparative biology.

Truncated radial glia as a common precursor in the late corticogenesis of gyrencephalic mammals.

The diversity of neural stem cells is a hallmark of the cerebral cortex development in gyrencephalic mammals, such as Primates and Carnivora. Among them, ferrets are a good model for mechanistic studies. However, information on their neural progenitor cells (NPC), termed radial glia (RG), is limited. Here, we surveyed the temporal series of single-cell transcriptomes of progenitors regarding ferret corticogenesis and found a conserved diversity and temporal trajectory between human and ferret NPC, despite the large timescale difference. We found truncated RG (tRG) in ferret cortical development, a progenitor subtype previously described in humans. The combination of in silico and in vivo analyses identified that tRG differentiate into both ependymal and astrogenic cells. Via transcriptomic comparison, we predict that this is also the case in humans. Our findings suggest that tRG plays a role in the formation of adult ventricles, thereby providing the architectural bases for brain expansion.

Elasmobranch genome sequencing reveals evolutionary trends of vertebrate karyotype organization.

Genomic studies of vertebrate chromosome evolution have long been hindered by the scarcity of chromosome-scale DNA sequences of some key taxa. One of those limiting taxa has been the elasmobranchs (sharks and rays), which harbor species often with numerous chromosomes and enlarged genomes. Here, we report the chromosome-scale genome assembly for the zebra shark Stegostoma tigrinum, an endangered species that has a relatively small genome among sharks (3.71 Gb), as well as for the whale shark Rhincodon typus Our analysis using a male-female comparison identified an X Chromosome, the first genomically characterized shark sex chromosome. The X Chromosome harbors the Hox C cluster whose intact linkage has not been shown for an elasmobranch fish. The sequenced shark genomes show a gradualism of chromosome length with remarkable length-dependent characteristics-shorter chromosomes tend to have higher GC content, gene density, synonymous substitution rate, and simple tandem repeat content as well as smaller gene length and lower interspersed repeat content. We challenge the traditional binary classification of karyotypes as with and without so-called microchromosomes. Even without microchromosomes, the length-dependent characteristics persist widely in nonmammalian vertebrates. Our investigation of elasmobranch karyotypes underpins their unique characteristics and provides clues for understanding how vertebrate karyotypes accommodate intragenomic heterogeneity to realize a complex readout. It also paves the way to dissecting more genomes with variable sizes to be sequenced at high quality.

Autocrine TGF-ß-positive feedback in profibrotic AT2-lineage cells plays a crucial role in non-inflammatory lung fibrogenesis.

The molecular etiology of idiopathic pulmonary fibrosis (IPF) has been extensively investigated to identify new therapeutic targets. Although anti-inflammatory treatments are not effective for patients with IPF, damaged alveolar epithelial cells play a critical role in lung fibrogenesis. Here, we establish an organoid-based lung fibrosis model using mouse and human lung tissues to assess the direct communication between damaged alveolar type II (AT2)-lineage cells and lung fibroblasts by excluding immune cells. Using this in vitro model and mouse genetics, we demonstrate that bleomycin causes DNA damage and activates p53 signaling in AT2-lineage cells, leading to AT2-to-AT1 transition-like state with a senescence-associated secretory phenotype (SASP). Among SASP-related factors, TGF-ß plays an exclusive role in promoting lung fibroblast-to-myofibroblast differentiation. Moreover, the autocrine TGF-ß-positive feedback loop in AT2-lineage cells is a critical cellular system in non-inflammatory lung fibrogenesis. These findings provide insights into the mechanism of IPF and potential therapeutic targets.

Single-cell transcriptome atlas of Drosophila gastrula 2.0.

During development, positional information directs cells to specific fates, leading them to differentiate with their own transcriptomes and express specific behaviors and functions. However, the mechanisms underlying these processes in a genome-wide view remain ambiguous, partly because the single-cell transcriptomic data of early developing embryos containing accurate spatial and lineage information are still lacking. Here, we report a single-cell transcriptome atlas of Drosophila gastrulae, divided into 77 transcriptomically distinct clusters. We find that the expression profiles of plasma-membrane-related genes, but not those of transcription-factor genes, represent each germ layer, supporting the nonequivalent contribution of each transcription-factor mRNA level to effector gene expression profiles at the transcriptome level. We also reconstruct the spatial expression patterns of all genes at the single-cell stripe level as the smallest unit. This atlas is an important resource for the genome-wide understanding of the mechanisms by which genes cooperatively orchestrate Drosophila gastrulation.

Squalomix: shark and ray genome analysis consortium and its data sharing platform.

The taxon Elasmobranchii (sharks and rays) contains one of the long-established evolutionary lineages of vertebrates with a tantalizing collection of species occupying critical aquatic habitats. To overcome the current limitation in molecular resources, we launched the Squalomix Consortium in 2020 to promote a genome-wide array of molecular approaches, specifically targeting shark and ray species. Among the various bottlenecks in working with elasmobranchs are their elusiveness and low fecundity as well as the large and highly repetitive genomes. Their peculiar body fluid composition has also hindered the establishment of methods to perform routine cell culturing required for their karyotyping. In the Squalomix consortium, these obstacles are expected to be solved through a combination of in-house cytological techniques including karyotyping of cultured cells, chromatin preparation for Hi-C data acquisition, and high fidelity long-read sequencing. The resources and products obtained in this consortium, including genome and transcriptome sequences, a genome browser powered by JBrowse2 to visualize sequence alignments, and comprehensive matrices of gene expression profiles for selected species are accessible through https://github.com/Squalomix/info.

Molecular and morphological investigations on the renal mechanisms enabling euryhalinity of red stingray Hemitrygon akajei.

Most cartilaginous fishes live in seawater (SW), but a few exceptional elasmobranchs (sharks and rays) are euryhaline and can acclimate to freshwater (FW) environments. The plasma of elasmobranchs is high in NaCl and urea concentrations, which constrains osmotic water loss. However, these euryhaline elasmobranchs maintain high levels of plasma NaCl and urea even when acclimating to low salinity, resulting in a strong osmotic gradient from external environment to body fluid. The kidney consequently produces a large volume of dilute urine to cope with the water influx. In the present study, we investigated the molecular mechanisms of dilute urine production in the kidney of Japanese red stingray, Hemitrygon akajei, transferred from SW to low-salinity environments. We showed that red stingray maintained high plasma NaCl and urea levels by reabsorbing more osmolytes in the kidney when transferred to low salinity. RNA-seq and qPCR analyses were conducted to identify genes involved in NaCl and urea reabsorption under the low-salinity conditions, and the upregulated gene expressions of Na+-K+-Cl- cotransporter 2 (nkcc2) and Na+/K+-ATPase (nka) were found in the FW-acclimated individuals. These upregulations occurred in the early distal tubule (EDT) in the bundle zone of the kidney, which coils around the proximal and collecting tubules to form the highly convoluted structure of batoid nephron. Considering the previously proposed model for urea reabsorption, the upregulation of nkcc2 and nka not only causes the reabsorption of NaCl in the EDT, but potentially also supports enhanced urea reabsorption and eventually the production of dilute urine in FW-acclimated individuals. We propose advantageous characteristics of the batoid-type nephron that facilitate acclimation to a wide range of salinities, which might have allowed the batoids to expand their habitats.

Repeated translocation of a supergene underlying rapid sex chromosome turnover in Takifugu pufferfish.

Recent studies have revealed a surprising diversity of sex chromosomes in vertebrates. However, the detailed mechanism of their turnover is still elusive. To understand this process, it is necessary to compare closely related species in terms of sex-determining genes and the chromosomes harboring them. Here, we explored the genus Takifugu, in which one strong candidate sex-determining gene, Amhr2, has been identified. To trace the processes involved in transitions in the sex-determination system in this genus, we studied 12 species and found that while the Amhr2 locus likely determines sex in the majority of Takifugu species, three species have acquired sex-determining loci at different chromosomal locations. Nevertheless, the generation of genome assemblies for the three species revealed that they share a portion of the male-specific supergene that contains a candidate sex-determining gene, GsdfY, along with genes that potentially play a role in male fitness. The shared supergene spans ∼100 kb and is flanked by two duplicated regions characterized by CACTA transposable elements. These results suggest that the shared supergene has taken over the role of sex-determining locus from Amhr2 in lineages leading to the three species, and repeated translocations of the supergene underlie the turnover of sex chromosomes in these lineages. These findings highlight the underestimated role of a mobile supergene in the turnover of sex chromosomes in vertebrates.

Evolution of variable lymphocyte receptor B antibody loci in jawless vertebrates.

Three types of variable lymphocyte receptor (VLR) genes, VLRA, VLRB, and VLRC, encode antigen recognition receptors in the extant jawless vertebrates, lampreys and hagfish. The somatically diversified repertoires of these VLRs are generated by serial stepwise copying of leucine-rich repeat (LRR) sequences into an incomplete germline VLR gene. Lymphocytes that express VLRA or VLRC are T cell-like, while VLRB-expressing cells are B cell-like. Here, we analyze the composition of the VLRB locus in different jawless vertebrates to elucidate its configuration and evolutionary modification. The incomplete germline VLRB genes of two hagfish species contain short noncoding intervening sequences, whereas germline VLRB genes in six representative lamprey species have much longer intervening sequences that exhibit notable genomic variation. Genomic clusters of potential LRR cassette donors, fragments of which are copied to complete VLRB gene assembly, are identified in Japanese lamprey and sea lamprey. In the sea lamprey, 428 LRR cassettes are located in five clusters spread over a total of 1.7 Mbp of chromosomal DNA. Preferential usage of the different donor cassettes for VLRB assemblage is characterized in our analysis, which reveals evolutionary modifications of the lamprey VLRB genes, elucidates the organization of the complex VLRB locus, and provides a comprehensive catalog of donor VLRB cassettes in sea lamprey and Japanese lamprey.

Technical considerations in Hi-C scaffolding and evaluation of chromosome-scale genome assemblies.

The recent development of ecological studies has been fueled by the introduction of massive information based on chromosome-scale genome sequences, even for species for which genetic linkage is not accessible. This was enabled mainly by the application of Hi-C, a method for genome-wide chromosome conformation capture that was originally developed for investigating the long-range interaction of chromatins. Performing genomic scaffolding using Hi-C data is highly resource-demanding and employs elaborate laboratory steps for sample preparation. It starts with building a primary genome sequence assembly as an input, which is followed by computation for genome scaffolding using Hi-C data, requiring careful validation. This article presents technical considerations for obtaining optimal Hi-C scaffolding results and provides a test case of its application to a reptile species, the Madagascar ground gecko (Paroedura picta). Among the metrics that are frequently used for evaluating scaffolding results, we investigate the validity of the completeness assessment of chromosome-scale genome assemblies using single-copy reference orthologues.

Chromosomal-scale de novo genome assemblies of Cynomolgus Macaque and Common Marmoset.

Cynomolgus macaque (Macaca fascicularis) and common marmoset (Callithrix jacchus) have been widely used in human biomedical research. Long-standing primate genome assemblies used the human genome as a reference for ordering and orienting the assembled fragments into chromosomes. Here we performed de novo genome assembly of these two species without any human genome-based bias observed in the genome assemblies released earlier. We assembled PacBio long reads, and the resultant contigs were scaffolded with Hi-C data, which were further refined based on Hi-C contact maps and alternate de novo assemblies. The assemblies achieved scaffold N50 lengths of 149 Mb and 137 Mb for cynomolgus macaque and common marmoset, respectively. The high fidelity of our assembly is also ascertained by BAC-end concordance in common marmoset. Our assembly of cynomolgus macaque outperformed all the available assemblies of this species in terms of contiguity. The chromosome-scale genome assemblies produced in this study are valuable resources for non-human primate models and provide an important baseline in human biomedical research.

Cooperation of LIM domain-binding 2 (LDB2) with EGR in the pathogenesis of schizophrenia.

Genomic defects with large effect size can help elucidate unknown pathologic architecture of mental disorders. We previously reported on a patient with schizophrenia and a balanced translocation between chromosomes 4 and 13 and found that the breakpoint within chromosome 4 is located near the LDB2 gene. We show here that Ldb2 knockout (KO) mice displayed multiple deficits relevant to mental disorders. In particular, Ldb2 KO mice exhibited deficits in the fear-conditioning paradigm. Analysis of the amygdala suggested that dysregulation of synaptic activities controlled by the immediate early gene Arc is involved in the phenotypes. We show that LDB2 forms protein complexes with known transcription factors. Consistently, ChIP-seq analyses indicated that LDB2 binds to > 10,000 genomic sites in human neurospheres. We found that many of those sites, including the promoter region of ARC, are occupied by EGR transcription factors. Our previous study showed an association of the EGR family genes with schizophrenia. Collectively, the findings suggest that dysregulation in the gene expression controlled by the LDB2-EGR axis underlies a pathogenesis of subset of mental disorders.

Developmental hourglass and heterochronic shifts in fin and limb development.

How genetic changes are linked to morphological novelties and developmental constraints remains elusive. Here, we investigate genetic apparatuses that distinguish fish fins from tetrapod limbs by analyzing transcriptomes and open-chromatin regions (OCRs). Specifically, we compared mouse forelimb buds with the pectoral fin buds of an elasmobranch, the brown-banded bamboo shark (Chiloscyllium punctatum). A transcriptomic comparison with an accurate orthology map revealed both a mass heterochrony and hourglass-shaped conservation of gene expression between fins and limbs. Furthermore, open-chromatin analysis suggested that access to conserved regulatory sequences is transiently increased during mid-stage limb development. During this stage, stage-specific and tissue-specific OCRs were also enriched. Together, early and late stages of fin/limb development are more permissive to mutations than middle stages, which may have contributed to major morphological changes during the fin-to-limb evolution. We hypothesize that the middle stages are constrained by regulatory complexity that results from dynamic and tissue-specific transcriptional controls.

Animals come in all shapes and sizes. This diversity arose through genetic mutations during evolution, but it is unclear exactly how these variations led to the formation of new shapes. There is increasing evidence to suggest that not all shapes are possible and that variability between animals is limited by a phenomenon known as "developmental constraint". These limitations direct parts of the body towards a specific shape as they develop in the embryo. Therefore, understanding the mechanisms underlying these developmental constraints could help explain how different body shapes evolved. The limbs of humans and other mammals evolved from the fins of fish, and this transition is often used to study the role developmental constraints play in evolution. This is an ideal model as there is already a detailed fossil record mapping this evolutionary event, and data pinpointing some of the genes involved in the development of limbs and fins. But this data is incomplete, and a full comparison between the genes activated in the fin and the limb during embryonic development had not been achieved. This is because most fish used for research have undergone recent genetic changes, making it hard to spot which genetic differences are linked to the evolution of the limb. To overcome this barrier, Onimaru et al. compared genetic data from the developing limbs of mice to the developing fins of the brown-banded bamboo shark, which evolves much slower than other fish. This revealed that although many genes commonly played a role in the development of the fin and the limb in the embryo, the activity of these shared genes was not the same. For example, genes that switched on in the late stages of limb development, switched off in the late stages of fin development. But in the middle of development, those differences were relatively small and both species activated very similar sets of genes. Many of these genes were pleiotropic, which means they have important roles in other tissues and therefore mutate less often. This suggests that the mid-stage of limb development is under the strongest level of constraint. Darwin's theory of natural selection explains that mutations drive evolution. But the theory cannot predict what kinds of new body shapes new mutations will produce. Understanding how the activity levels of different genes affect development could help to fill this knowledge gap. This has potential medical applications, for example, understanding why some genetic changes cause more serious problems than others. This work suggests that mutations in genes that are active during the mid-stage of limb development may have the most serious impact.

Facilitated NaCl Uptake in the Highly Developed Bundle of the Nephron in Japanese Red Stingray Hemitrygon akajei Revealed by Comparative Anatomy and Molecular Mapping.

Batoidea (rays and skates) is a monophyletic subgroup of elasmobranchs that diverged from the common ancestor with Selachii (sharks) about 270 Mya. A larger number of batoids can adapt to low-salinity environments, in contrast to sharks, which are mostly stenohaline marine species. Among osmoregulatory organs of elasmobranchs, the kidney is known to be dedicated to urea retention in ureosmotic cartilaginous fishes. However, we know little regarding urea reabsorbing mechanisms in the kidney of batoids. Here, we performed physiological and histological investigations on the nephrons in the red stingray (Hemitrygon akajei) and two shark species. We found that the urine/plasma ratios of salt and urea concentrations in the stingray are significantly lower than those in cloudy catshark (Scyliorhinus torazame) under natural seawater, indicating that the kidney of stingray more strongly reabsorbs these osmolytes. By comparing the three-dimensional images of nephrons between stingray and banded houndshark (Triakis scyllium), we showed that the tubular bundle of stingray has a more compact configuration. In the compact tubular bundle of stingray kidney, the distal diluting tubule was highly developed and frequently coiled around the proximal and collecting tubules. Furthermore, co-expression of NKAα1 (Na+/K +-ATPase) and NKCC2 (Na+- K+-2Cl- cotransporter 2) mRNAs was prominent in the coiled diluting segment. These findings imply that NaCl reabsorption is greatly facilitated in the stingray kidney, resulting in a higher reabsorption rate of urea. Lowering the loss of osmolytes in the glomerular filtrate is likely favorable to the adaptability of batoids to a wide range of environmental salinity.

Novel developmental bases for the evolution of hypobranchial muscles in vertebrates.

BACKGROUND: Vertebrates are characterized by possession of hypobranchial muscles (HBMs). Cyclostomes, or modern jawless vertebrates, possess a rudimentary and superficial HBM lateral to the pharynx, whereas the HBM in jawed vertebrates is internalized and anteroposteriorly specified. Precursor cells of the HBM, marked by expression of Lbx1, originate from somites and undergo extensive migration before becoming innervated by the hypoglossal nerve. How the complex form of HBM arose in evolution is relevant to the establishment of the vertebrate body plan, but despite having long been assumed to be similar to that of limb muscles, modification of developmental mechanisms of HBM remains enigmatic. RESULTS: Here we characterize the expression of Lbx genes in lamprey and hagfish (cyclostomes) and catshark (gnathostome; jawed vertebrates). We show that the expression patterns of the single cyclostome Lbx homologue, Lbx-A, do not resemble the somitic expression of mammalian Lbx1. Disruption of Lbx-A revealed that LjLbx-A is required for the formation of both HBM and body wall muscles, likely due to the insufficient extension of precursor cells rather than to hindered muscle differentiation. Both homologues of Lbx in the catshark were expressed in the somitic muscle primordia, unlike in amniotes. During catshark embryogenesis, Lbx2 is expressed in the caudal HBM as well as in the abdominal rectus muscle, similar to lamprey Lbx-A, whereas Lbx1 marks the rostral HBM and pectoral fin muscle. CONCLUSIONS: We conclude that the vertebrate HBM primarily emerged as a specialized somatic muscle to cover the pharynx, and the anterior internalized HBM of the gnathostomes is likely a novelty added rostral to the cyclostome-like HBM, for which duplication and functionalization of Lbx genes would have been a prerequisite.

Early vertebrate origin of CTCFL, a CTCF paralog, revealed by proximity-guided shark genome scaffolding.

The nuclear protein CCCTC-binding factor (CTCF) contributes as an insulator to chromatin organization in diverse animals. The gene encoding this protein has a paralog which was first identified to be expressed exclusively in the testis in mammals and designated as CTCFL (also called BORIS). CTCFL orthologs were reported only among amniotes, and thus CTCFL was once thought to have arisen in the amniote lineage. In this study, we identified elasmobranch CTCFL orthologs, and investigated its origin with the aid of a shark genome assembly improved by proximity-guided scaffolding. Our analysis employing evolutionary interpretation of syntenic gene location suggested an earlier timing of the gene duplication between CTCF and CTCFL than previously thought, that is, around the common ancestor of extant vertebrates. Also, our transcriptomic sequencing revealed a biased expression of the catshark CTCFL in the testis, suggesting the origin of the tissue-specific localization in mammals more than 400 million years ago. To understand the historical process of the functional consolidation of the long-standing chromatin regulator CTCF, its additional paralogs remaining in some of the descendant lineages for spatially restricted transcript distribution should be taken into consideration.

TGF-ß1-driven reduction of cytoglobin leads to oxidative DNA damage in stellate cells during non-alcoholic steatohepatitis.

BACKGROUND & AIMS: Cytoglobin (CYGB) is a respiratory protein that acts as a scavenger of reactive oxygen species. The molecular role of CYGB in human hepatic stellate cell (HSC) activation and human liver disease remains uncharacterised. The aim of this study was to reveal the mechanism by which the TGF-ß1/SMAD2 pathway regulates the human CYGB promoter and the pathophysiological function of CYGB in human non-alcoholic steatohepatitis (NASH). METHODS: Immunohistochemical staining was performed using human NASH biopsy specimens. Molecular and biochemical analyses were performed by western blotting, quantitative PCR, and luciferase and immunoprecipitation assays. Hydroxyl radicals (•OH) and oxidative DNA damage were measured using an •OH-detectable probe and 8-hydroxy-2'-deoxyguanosine (8-OHdG) ELISA. RESULTS: In culture, TGF-ß1-pretreated human HSCs exhibited lower CYGB levels - together with increased NADPH oxidase 4 (NOX4) expression - and were primed for H2O2-triggered •OH production and 8-OHdG generation; overexpression of human CYGB in human HSCs reversed these effects. Electron spin resonance demonstrated the direct •OH scavenging activity of recombinant human CYGB. Mechanistically, pSMAD2 reduced CYGB transcription by recruiting the M1 repressor isoform of SP3 to the human CYGB promoter at nucleotide positions +2-+13 from the transcription start site. The same repression did not occur on the mouse Cygb promoter. TGF-ß1/SMAD3 mediated αSMA and collagen expression. Consistent with observations in cultured human HSCs, CYGB expression was negligible, but 8-OHdG was abundant, in activated αSMA+pSMAD2+- and αSMA+NOX4+-positive hepatic stellate cells from patients with NASH and advanced fibrosis. CONCLUSIONS: Downregulation of CYGB by the TGF-ß1/pSMAD2/SP3-M1 pathway brings about •OH-dependent oxidative DNA damage in activated hepatic stellate cells from patients with NASH. LAY SUMMARY: Cytoglobin (CYGB) is a respiratory protein that acts as a scavenger of reactive oxygen species and protects cells from oxidative DNA damage. Herein, we show that the cytokine TGF-ß1 downregulates human CYGB expression. This leads to oxidative DNA damage in activated hepatic stellate cells. Our findings provide new insights into the relationship between CYGB expression and the pathophysiology of fibrosis in patients with non-alcoholic steatohepatitis.

Multifaceted Hi-C benchmarking: what makes a difference in chromosome-scale genome scaffolding?

BACKGROUND: Hi-C is derived from chromosome conformation capture (3C) and targets chromatin contacts on a genomic scale. This method has also been used frequently in scaffolding nucleotide sequences obtained by de novo genome sequencing and assembly, in which the number of resultant sequences rarely converges to the chromosome number. Despite its prevalent use, the sample preparation methods for Hi-C have not been intensively discussed, especially from the standpoint of genome scaffolding. RESULTS: To gain insight into the best practice of Hi-C scaffolding, we performed a multifaceted methodological comparison using vertebrate samples and optimized various factors during sample preparation, sequencing, and computation. As a result, we identified several key factors that helped improve Hi-C scaffolding, including the choice and preparation of tissues, library preparation conditions, the choice of restriction enzyme(s), and the choice of scaffolding program and its usage. CONCLUSIONS: This study provides the first comparison of multiple sample preparation kits/protocols and computational programs for Hi-C scaffolding by an academic third party. We introduce a customized protocol designated "inexpensive and controllable Hi-C (iconHi-C) protocol," which incorporates the optimal conditions identified in this study, and demonstrate this technique on chromosome-scale genome sequences of the Chinese softshell turtle Pelodiscus sinensis.

TChIP-Seq: Cell-Type-Specific Epigenome Profiling.

Epigenetic regulation plays central roles in gene expression. Since histone modification was discovered in the 1960s, its physiological and pathological functions have been extensively studied. Indeed, the advent of next-generation deep sequencing and chromatin immunoprecipitation (ChIP) via specific histone modification antibodies has revolutionized our view of epigenetic regulation across the genome. Conversely, tissues typically consist of diverse cell types, and their complex mixture poses analytic challenges to investigating the epigenome in a particular cell type. To address the cell type-specific chromatin state in a genome-wide manner, we recently developed tandem chromatin immunoprecipitation sequencing (tChIP-Seq), which is based on the selective purification of chromatin by tagged core histone proteins from cell types of interest, followed by ChIP-Seq. The goal of this protocol is the introduction of best practices of tChIP-Seq. This technique provides a versatile tool for tissue-specific epigenome investigation in diverse histone modifications and model organisms.